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A system consisting of two identical single-mode cavities coupled to a common environment is 
investigated within the framework of algebraic dynamics. Based on the left and right representations 
of the Heisenberg-Weyl algebra, the algebraic structure of the master equation is explored and exact 
analytical solutions of this system are obtained. It is shown that for such a system, the environment 
can produce entanglement in contrast to its commonly believed role of destroying entanglement. In 
addition, the collective zero-mode eigen solutions of the system are found to be free of decoherence 
against the dissipation of the environment. These decoherence-free states may be useful in quantum 
information and quantum computation. 
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I. INTRODUCTION 



Quantum entanglement plays a basic role in quantum communication and quantum computation [1]. The creation 
of entanglement between qubits is of fundamental importance for further quantum computation processing. The 
\ entanglement can be created by a direct interactions between qubits [2] or an indirect interactions via a third party [3] . 
However, both of the above processes are confined in closed system, i.e., the influences of environment are neglected. 
Recent investigations showed that environment can be helpful to the entanglement creation in an open system [4, 5], 
which provide a perspective to use the environment to implement decoherence-free quantum information processing 
[6, 7]. In order to treat the influences of the environment on an open quantum system, the Born-Markovian master 
, equation approach has been widely used. The common feature of the quantum master equations is the existence of the 
sandwich terms of the Liouville operator where the reduced density matrix of the system is in between some quantum 
excitation and de-excitation operators. These terms result from the elimination of an enormous irrelevant degrees 
of freedom of the environment. Except for some simple cases, for example, a single-mode of the cavity field coupled 
to the vacuum or stationary regime properties [8], it is very difficult to solve directly the master equation. Instead, 
it is usual to convert the master equations into some c-number equations in the coherent state representation-the 
Fokker-Planck equation [9, 10]. 

In the previous works [11-13], we have proposed and developed an algebraic method to treat the sandwich terms in 
the Liouville operator for quantum statistical systems. This method is just a generalization of the algebraic dynamical 
method [14] from quantum mechanical systems to quantum statistical systems. According to the characteristic of 
the sandwich terms in the Liouville operator, the left and right representations of the relevant algebra [15] have been 
introduced and the corresponding composite algebra has been constructed. As a result, the master equation has been 
converted into a Schrodinger-like equation and the problems can be solved exactly. 

In this paper, we shall use this method to solve the problem of two identical cavities coupled to a common en- 
vironment. The system consisting of two coupled cavities and a similar system of two coupled harmonic oscillators 
are important in quantum optics and quantum information theory. From these systems, many properties, such as, 
information transfer of quantum states [16], and the quantum statistical properties of the two-coupled modes of elec- 
tromagnetic fields [17], have been investigated. In Ref. [18], the entanglement of the two coupled harmonic oscillator 
system has been quantitatively studied. However, all these studies neglect the effect of environment on the modes 
of cavities and the system is thus closed. Although Ref. [19] studied the dissipative two oscillators system, it was 
based on the quantum characteristic function approach and the direct analytical solutions are not obtained. In this 
work, we shall consider a specific coupling of the two cavity fields, induced through the individual interactions of 
the two cavities with the environment and our system is thus different from Ref. [19]. By introducing a collective 
mode consisting of the two cavity modes, the quantum master equation is obtained. The two independent cavities 
are thus interacting with each other indirectly through the environment. Using the algebraic dynamical method, the 
full algebraic structure of the master equation is explored and its exact analytical solution is obtained. From the zero 
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eigenvalue solutions, the decoherence-free sub-space is obtained. Finally, the exact entanglement dynamics for the 
single- and two-photon processes is obtained from the time-dependent solutions. We find that the system exhibits an 
interesting feature of environment- assisted entanglement generation. Our present study shows that the environment 
is not just a negative source of decoherence, it also can be a positive source of entanglement. 

The paper is organized as follows. In Sec. II, the model Hamiltonian of the system is presented and the master 
equation for the reduced density matrix of the system is obtained. In Sec. Ill, based on the left and right represen- 
tations of the relevant algebra, we introduce a composite algebra, in terms of which the dynamical su(l, 1) plus u{2) 
algebraic structures of the Liouville operator( rate operator ) of the master equation are found. Sec. IV is devoted to 
the analytical solutions of the master equation. The applications of its time-dependent solutions to the entanglement 
dynamics for the single- and two-photon processes are studied in Sec. V. Finally, a brief summary is given in Sec. VI. 



II. TWO IDENTICAL CAVITIES IN A COMMON ENVIRONMENT 



Consider two identical cavities (or two harmonic oscillators) interacting with a common environment. With the 
dipole interaction and in the rotating wave approximation, the system can be described by the Hamiltonian 



H = lo a \ a i + ^kblbk + y^(ffifcqtfrfc + H.c.) 



i k i,k 

where di(i = 1,2) arc boson operators of the two cavities with the same frequency lo and bk are those for the 
environment fields, respectively, gik is the coupling constant. It is reasonable to assume that the coupling of the 
two cavities to the environment is same, i.e., gik = gk- This is a generation of Dicke limit of quantum optics [20]. 
Using the standard technique of quantum optics, one can obtain the master equation for the reduced density of the 
two-cavity system under the standard Born-Markovian approximation [8, 10] 



dp(t) 



5-{[2aip(t)a\ - a\aip{t) - p{t)a\a{\ + [2a 2 p{t)a\ - a\a 2 p(t) - p{t)a\a 2 ]} 



dt 2 

+^{[2a\p{t)a\ - a[a 2 p(t) - p{t)a\a 2 ] + [2a 2 p{t)a\ - a\a\p{t) - p{t)a\a{\}, (1) 

where q{> 0) is the decay constant of the collective mode. It is noted that on the right-hand side of Eq. (1) the first two 
terms denote the individual dissipation of the two cavities due to the environment, while the last two terms describe 
the coupling (photon exchange) between the two cavities mediated by the environment. Eq. (1) is complicated and 
can be simplified by introducing a collective (quasi-photon) mode operator A = -^=(ai + a 2 ), as did in Ref. [21]. In 
doing so, the master equation of the reduced density becomes 

^ = \[2Ap{t)rf dAp® p(t)A*A]. (2) 
In the following we shall discuss the algebraic structure of this system. 



III. ALGEBRAIC STRUCTURE OF THE MASTER EQUATION 



Based on the left and right representations of Heisenberg-Weyl algebra [12], the master equation (2) can be converted 
into a Schrodinger-like equation and the algebraic structure can thus be established. Noticing the fact that the collec- 
tive operators A and A^ obey the same commutation relations as the photon operators a and a) do, {A, A\ A^ A, 1} 
thus constitutes a /iu>(4) algebra. Define the left and right algebras of /in; (4) 

hw{A) r = {A\A r \N r = A rt A r ,l}, 

hw(4)i = {A\A l \N l = A l ^A\l}, (3) 

where hw(4) r acts to the right on the ket photon number state \n) and hw(4)i acts to the left on the bra photon 
number state (n\. They have the following commutators 



[A r ,A rJ< } = l,[A r ,N r ] = A r ,[A r \N r ] = -A rt ; 
[A l ,A lJf ] = -l,[A l ,N l ] = -A l ,[A l ^,N l ] = A l l 



(4) 
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It is noted that hw(A)i (hw(4) r ) is isomorphic (anti-isomorphic) to hw(4). Since hw(4)i and hw(A)i act on different 
spaces-bra and ket spaces, the operators commute with each other, i.e., [hw(4)i, hw(A) r ] = 0. The reduced density 
matrix is a vector of the von Nuemann super space, which has the hw(4)i ® hw(4) r algebraic structure containing 
sp(4) as a relevant sub-algebra, as shown below. 

From the above basic algebras we can constitute a composite algebra C in the adjoint representation, 

C = {K_ = A r A l \K + = A r ^A l ,K = ^(N r + N 1 )} 

We see that C is an su(l, 1) algebra satisfying the following commutation rules 

[K ,K ± ] = ±K ± ,[K_,K + ] = 2K , 

which is an sub-algebra of the sp{4) and can be derived from Eqs. (4). The operator K_ acts on the photon number 
bases of the von Neumann space as follows 

K-\nin 2 ){mim 2 \ = ^(\/m|m - 1, n 2 ) + \fm\n\, n 2 - 1)) 

((mi - l,m 2 \*Jm~[ + (mi,m 2 - ll^/m^). (5) 

Similarly, for the super- vector bases of the quasi-photon (collective photon) number states \n) qq (m\ marked by the 
subscript q, the actions of the de-excitation operator K- and the number operator K n read 

K-\n) qq (m\ — y/nm\n — l) qq (m — 1|, 
is i \ / i n + m + 1 

K \n} qq (m\ = \n) qq {m\. (6) 

With this composite algebra at hand, it is straightforward to convert the master equation (2) into Schrodinger-like 
equation 

d ~f - ■>«>• 

r = ?Jf_ - <;K + i (7) 

Since the rate operator T is a linear function of the su(l, 1) generators, we conclude that the master equation possesses 
an su(l, 1) dynamical symmetry in the quasi-photon number representation. Thus the system is integrable and can 
be solved analytically according to the algebraic dynamics [14]. 

It is remarkable to note that Kq is diagonalized in the quasi-photon representation, but it is not so in the real 
photon number representation of the von Neumann space of the density matrix. As will be seen soon, it has a 
different algebraic structure in the real photon number representation. In terms of the two real photon operators, K 
can be written as 

K = ^(N 1 +N 2 + S + + S_), 

where Ni — a\ + d} a\ = n \ + n\{i = 1, 2), S + = a 2 + a^a^, and S- = a^a^ + a l ^a l 2 . It is not difficult to find 
that the operators (Ni,N 2 , S+, S-) form a u(2)(— u(l) + su(2)) algebra 

u(2) = {.So = (m - N 2 )/2,S+,S-;N = (N 1 +N 2 )/2}, 

which obey the following commutation rules 

[N, S±] - [N, S ] = 0, 

[S ,S±] = ±S ± ,[S-,S + } = -2So. 

These operators act on the von Neumann space in the real photon number representation as follows: 



S + \nin 2 )(mim 2 \ — \J (m + l)n 2 |ni + 1, n 2 — l)(mim 2 \ + \J (mi + \)m 2 \n\n 2 ) (mi + 1, m 2 — 1|, 

5_|n 1 n 2 )(TOim 2 | = \/ni(n 2 + l)|ni - l,n 2 + 1)(toito 2 | + y / m 1 (m 2 + l)\n 1 n 2 )(mi - l,m 2 + 1|, 

Ni\nin 2 )(mim 2 \ = (n 4 + rrn + l)\nin 2 )(m\m 2 \, (i = 1, 2) 

So\nin 2 ) (m\m 2 \ = (ni + mi - n 2 - m 2 )/2\nin 2 )(mim 2 \, 

N\mn 2 )(m\m 2 \ = (ni + mi + n 2 + m 2 + 2)/2|nin 2 )(mim 2 |. (8) 
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Up to now we have explored the complete dynamical symmetry of the system and the results can be summarized 
as follows: the largest dynamical algebra of the two-cavity system is the hw(4)i (g> hw(4) r D sp(4) algebra since it 
has two kinds of photons; due to the special structure of the rate operator T, the system in fact has the su(l, 1) and 
u{2) sub-algebraic dynamical symmetries of the largest dynamical symmetry algebra sp(4). For the quasi-photon, the 
dynamical symmetry is of the su(l, 1) sub-algebra and for the two kinds of real photons it is u{2) algebra. 



IV. SOLUTIONS OF THE MASTER EQUATION 

A. Eigen solution and decoherence-free states 

To better understand the property of steady solution of the master equation, we first investigate the eigensolutions 
of the master equation. The eigen equation reads 

Tp = 1P . (9) 

To solve Eq. (9), we introduce two similarity transformations corresponding to the above two sub-algebraic structures. 
The first one is 

p = U lPl ,U 1 =e- K -. 

The operator T is thus diagonalized in the quasi-photon number representation as follows 

T 1 = Ur 1 rUi = -<;K + l (10) 
From Eq. (10) the zero-mode eigensolution of T in the quasi-photon number representation can be obtained 

po = |omo|, (ii) 

which has the same form as the stationary solution of the master equation for the single-mode damped harmonic 
oscillator [8]. 

One notes that in the quasi-photon number representation, the collective quasi-photon mode is used, in which two 
cavities are coupled. In experiments, it is convenient to use the real photon number representation. Moreover, in this 
representation, the interaction of the two cavities can also be studied. 

To diagonalize the rate operator f i(or K ) in the real photon number representation, the second similarity trans- 
formation is needed, namely, 



where U 2 is defined as 



Under the conditions 



f 2 = E/^IWa, (12) 
Pi = U 2P2 ,U 2 = e^ s +e^ s -. 



1-/3+ = 0, 

2/3+/?_ + l = 0, (13) 

i.e., /?_)_ = ±1 and /?_ = the rate operator f 2 is diagonalized in terms of both of the two Cartan operators N and 
S of the u(2) algebra and the real photon number representation^!, N 2 ), 

= ^(f2N 1+ f 3 N 2 ) + ^ 

h = 

From the expression of T 2 we can see that the eigen solutions of T 2 are highly degenerate with respect to different 
photon number distributions among the bra and ket states. So besides the two Cartan operators we need two additional 
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quantum numbers to distinguish the degenerate states. Then the complete set of commutation quantum operators 
is {iV, So, n^, n 2 }. For convenience, we choose the equivalent set of the complete commutation quantum operators as 
{n\, n{,n r 2 ,n l 2 }. 

After making two inverse transformations of T 2 and T\, we obtain the eigen solution of Eq. (9) on the real photon 
number bases of the von Neumann space as follows, 

7(nin 2; mim 2 ) = -^[/ 2 (rti + mi + 1) + / 3 (n 2 + m 2 + 1)] + |, 

p(n 1 n 2 .m 1 m 2 ) = e~ K ~ e^ +s+ e^~ s - |nin 2 )(mira 2 |. (14) 

Substituting the two sets of solutions of Eqs. (13) into the expression of 7, we see that both of them contain the 
quasi-photon zero-mode solution. For the second similarity transformation, the two sets of solutions of Eqs. (13) 
are physically equivalent and degenerate. This can be seen from the expressions of 7, / 2 , and fa in terms of the 
transformation coefficient /3+ which produces the solutions labeled by n\ and m\ (or n 2 and rn 2 ), since / 2 and f$ 
cannot be non-zero simultaneously under the similarity transformation. The degeneracy of the zero quasi-photon 
mode comes from the different real photon number distributions among the two cavities. From Eqs. (14) we can 
express the zero quasi-photon mode solution of the system in terms of the real photon number bases as follows 

po = c e- K -e s +e'^ s -\nO)(mO\ = c' e- K - e- s +e^ s - \0n)(0m\ (15) 

where c and c' are normalized constants. After some straightforward calculations, it is found that the action of 
er K - have no effect on the states e s+ e~z s - |nO)(mO| (or e~ s +e2 S - |0n)(0m|). Since Tp = 0, po arc an invariant 
subspace under the time evolution of the master equation (7). It is found that the quasi-photon zero- mode subspace 
(15) contains many entangled real photon states in the two cavities, which are thus decohence-free with respect to the 
action of the dissipative operator T in the master equation (7). For example, for the zero-photon process, the zero-mode 

state is p °^ = |00)(00|; for the single-photon process, the zero-mode state is p^ = \<f>)(<f>\, where \<j>) = -^(|10) — 1 01) ) ; 

for the two-photon process, the zero-mode state is p 2 ^ = \tp)(ip\, where \ip) = |(|02) + 1 20) — \/2|ll)), etc. All the 
above entangled real photon states are the members of the quasi-photon zero-mode subspace and decoherence-free 
from the dissipative action of T in the master equation(7). 

From the above analysis, we see that the quasi-photon zero-mode solution (15) forms a highly degenerate and 
dissipation-free collective subspace in terms of different real photon number states which are orthogonal to each other, 
and many of them are entangled and thus decoherence-free. It is noted that the equilibrium steady solutions of the 
system depending on initial states are not unique, as shown later. All of them are within the quasi-photon zero-mode 
subspace and consisting of all possible mixture of these real photon states with collective zero-mode. 



B. Time-dependent solutions of the master equation 



Next we shall investigate the dynamical properties of the system based on the real photon number representation. 
First we solve the master equation (7) and get its time-dependent solutions. Then we analyse the time evolution 
behavior of the solutions and prove that for some initial product states the entanglement can be produced by the 
environment induced dynamics and for some initial state the system is decoherence-free. 

By introducing a time-dependent gauge transformation 

P = U 1 (t)p, U 1 (t) = e a -W K -, 
the master equation can be rewritten as a diagonal form in the quasi-photon number representation 

dffl 



dt 



f = U ^(t)TU 1 (t)-U 1 -\t)^ = -<;K + i, 



if the following gauge transformation conditions are satisfied 

da-(t) 



di =?[l + a_(i)]. (16) 
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The time-dependent solution is now 

p(t) = U 1 (t)p(t) = Ui(t)e-lo^ d - p (0) 

= Y, C ^e a - {t)K ~e =5i ^\n) qq {rn\, (17) 

m.n 

where we have used the initial super state as follows p(0) = p(0) = n Cm,n\n}q q (m\. 

Similar to the steady case, it is necessary to express the above solution in the real-photon number representation. 
To this end, one needs the second time-dependent gauge transformation 

U 2 (t) = e^ s +e^ s - 

where the initial condition is taken as 1/2(0) = 1. Under the diagonalization conditions 



dt 4 
d/3_(i) -? 



(l + 2/? + (t)/?_(*)), (18) 



dt 4 

the rate operator T'(t) (thus T(tj) can be diagonalized in the real photon representation and reads 

f = (t) f'U 2 (t) (t) = ^[N- P + (t)S ] + \. 

After the two inverse transformations of U\(t) and ^(i), the time-dependent solution of Eq. (1) is then 

P {t) = u.w^mt) 

= m(t)U 2 (t) ^ 1 ^,m 1 ^ 2 e / « { ^ [(1 ^ +(r))( " 1+mi+1)+(1+ ^ (T))( " 2+m2+1)1+i}d 1n 1 n 2 )(m 1 m 2 |, (19) 

mi ,711 
n 2 ,m 2 

where we have set the initial state of the system as p(0) = p(0) = J2 mi ni n 2 m 2 Cn 1 ,n 2 ,m 1 ,m 2 \n\n2){ra\ra2\ in the real 
photon number representation. This solution is our central result in the present work. In the following we consider 
an explicit initial state to study explicitly the properties of this solution. 

V. APPLICATION OF THE SOLUTION 

In this section, we first consider explicitly the applications of the solution obtained above to the single- and two- 
photon processes, which show some interesting properties. Then we give a brief discussion on the implication to the 
quantum communication and quantum computation. 

A. Single-photon process 



Consider an initial state as p(0) = |x)(xl> where |%) = a|01) + \/l — a 2 |10), and a denotes the relevant amplitude. 
In terms of Eqs. (5, 8), the time-dependent solution (19) can be written as 

p(t) = Poi,oi(*)|01><01| + />io,io(*)|10><10| + /v/(u[|10)(01| + |01)(10|] + p O o,oo(U|00)(00|, (20) 

where the time-dependent coefficients can be obtained analytically, 

rt CO+(T) + l) , C(g+(T)-1) 

>0 2 dT , it 2\ n / 1 \ 2 10 



Poi,oi(*) = oV°" 2 T + (l-a 2 )(5_(t) 2 e° 2 " + 2a\/l^a? j3-(t)e 



t C(/3+(T) + l) t C(g+(T)-1) 

J o 2 " T . /-. 2^^1 1 a tx\o ^\i2 J o 



Pio.ioW = a 2 p + {tfe '°~ 2 T + (l-a 2 )[l + /3 + (t)/3_(t)] 2 e° 2 " + 2a^l^p + {t)[l + /3+(t)/3-(t)]e 

f t c(o+( T )+D j_ t c(g+(r)-i) _ s 

p o// (t) = a 2 f3+(t)e° 2 + (1 - a 2 )[l + (3 + (t)(3-(t)]p-(t)e° 2 + a^l - a 2 [l + 2/J+(t)/3_(t)]e 1 

Poo,oo(U = ^#{a 2 [l + /3 + (*)] 2 e° ^ +(l- a 2 )[l + /3_(t)+/3 + (U/3-(u] 2 e /0 ^^ 



+2av / l-« 2 [l + /?+ (*)] [1 + /?- (t) + /?+ (*)/?- (t)]e T }. (21) 
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From the above solution we can see that C/2 (i) makes the re-distribution of the photon number between the two 
cavities, which produces entanglement between them, while U\(t) induces de-excitation of the collective modes and 
leads to the zero collective mode. In order to study explicitly the time evolution of the system, it is needed to solve 
the time-dependent transformation parameters a_(t) and /3±(t), which can be obtained from Eq. (16) and Eqs. (18), 

a_(t) - e rt -l, 

/?+(*) = -tanh(|), 

/?-(*) = -^ sinh (f)- ( 22 ) 

It is known that the entanglement between the two subsystems can be measured by a quantity named as the entan- 
glement of formation [22]. This quantity can also be calculated by concurrence C(p) [23]. As usual, C{p) is defined 
as 

C(p) = max(0, Ai — A2 — A3 — A 4 ), 
where Aj(i = 1, • • • , 4) are eigenvalues of matrix 

[p 1 / 2 (a y ®a y )p*(a y ®a y )p 1 / 2 ] 1 / 2 

and Ai > A2 > A3 > A4. Here a y = ^ ^ ^ ^ and p* is the complex conjugation of p relative to the eigenbasis of <r z . 

It can be shown that the concurrence varies from C = for a disentangled state to C = 1 for a maximally entangled 
state. According to the solution (21) it is easy to find that the concurrence C(p) — 2\p ff(t)\ in the present case. 

In the following we will consider three cases of the entanglement and decoherence dynamics. For the first case 
we show that, for product state, the steady entanglement between the two subsystems can be produced by the 
environment. In the second case we present the decoherence-free state which is stable against the dissipation. The last 
case represents the situation where the environment plays the conventional role on system, i.e. it induces decoherence 
and makes the system disentangled. 

Case 1, a = 0. 

In this case the initial state is a product state and has no entanglement. Substituting Eqs. (22) into Eq. (21) we 
have 

Pit) = (1 ~ e 4 ~" )2 l01)(01| + (1 + e 4 ~" )2 |10)(10l + ^_l[|oi)<10| + |10)(01|] + IZ|Z!|00}<00|. 

After a long time-evolution, it is noted that all time-dependent coefficients become constant, and the density matrix 
can be denoted as 



\ 



Ps = 



\ x -I 



- 2 P + 2 P > 



\ 0/ 

p^ = |00)(00|, pW = \<t>)(<t>\ (23) 

which is in the basis of (|00), |01), |10), |11)} and pf(i = 0, 1) are obtained in section (IV. A) as zero-mode eigen 
solutions for zero- and single- photon processes, respectively. One can evaluate the concurrence of this steady state as 
C(p) = 0.5. The steady state is the possible mixture of the degenerate zero-mode eigen solutions. This is consistent 
with the discussion in section (IV. A). From the purification scheme proposed in Ref. [22] one can get the maximal 
entangled state \4>) from a collective pairs of the state above with probability 1/16. Thus it is shown that the maximal 
entangled state \<j>) between the two cavities can be produced from a product state by the environment. 

The time evolution behaviors of the density matrix and the concurrence are shown in Fig. 1 and 2, respectively. 

The same discussion can be applied to the initial product states with a — ±1, in which p(t) and C(p) have the 
similar time evolution behaviors. 

Case 2, a = -l/\/2 

The initial state is the maximally entangled state. In this case we found that the density matrix is time-independent, 
i.e., 

P (t) = i|01)(01| + ||10><10| - J[|01><10| + |10)(01|] = 
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Thus, this solution is decoherence-free and the concurrence is also a maximal constant C{p) = 1. From section (IV. 
A) we know that \4>){4>\ is the zero- mode eigen solution for single-photon process of the master equation (1). 
Case 3, a= l/y/2 

In this case we obtain the time dependent solution as follows 

P(t) = ^{|01)(01| + |10)(10| + [|01)(10| + |10)(01|]} + (1 - e-*)|00><00|. 

Its time evolution is plotted in Fig. 3, and asymptotically the state approach to a product state with zero real photon 
number. 



The concurrence C(p) decays quickly, as shown in Fig. 4, indicating that the decoherence is very strong and the two 
cavities will lose all of the information contained in its entangled initial state. 

Except for the above three special cases, it is difficult to obtain a simple expression of the density matrix for general 
cases of a and numerical calculations are needed. In Fig. 5, we show the time evolution of the density matrix for a 
general initial state characterized by a. Fig. 6 presents the time behavior of the corresponding concurrence. 



B. Two-photon process 

In this two-photon process, each mode of the two cavities is related to three states, |0), |1), and |2), i.e., a qutrit. 
Therefore, this system is a potential candidate for a qutrit quantum information processing [24, 25]. The advantage 
of qutrits instead of qubits is that the qutrits are much secure against the symmetric attacks in a quantum key 
distribution protocol, as shown by Bruss and Macchiavcllo [25]. In a recent experiment [26], the arbitrary qutrit 
states have been realized on the single-mode biphoton field, which shares some similarities with our model considered 
here. Based on the exact solution Eq. (19) in the two-photon process we shall analyse the decoherence effects on the 
qutrit state generation. 

Following Ref . [27] , the entanglement of bipartite qutrit states can be measured in terms of the Schmidt coefficients 
by 

1*3x3) = h\xi,yi) + k 2 \x 2 ,V2) + k 3 \x 3 ,y 3 ), 
<7(|$3x3» - \fmkl + k~M + k 2 k 2 3 ), (24) 

where Eq. (24) is a generalized definition of the concurrence for qutrit states. 

Suppose that the initial pure state is p(0) = , where \u>) — a\02) + bj 11) + c|20) where the coefficients a, 6, 

and c satisfy \J a 2 + b 2 + c 2 = 1. It is very difficult to discuss an arbitrary case. Instead, we focus on the following 
special cases. 

Case 1, a = l,b = c = 

The time-dependent solution can be obtained 

fe~ ? * + 1) 2 1 - e~ 2? * (1 - e" ? M 2 

P(t) = { 4 ' Mt))(<p(t)\ + — g — W)){<Kt)\ + - — i — ^-|oo><oo|, 

(l + e'"*/ 2 ) 2 1-e"^ (l-e'"^ 2 ) 2 

= Ww |02> - TSF^TI) 111 ' + Wt>°>' 

1 + e-^l 2 1 - e-^l 2 

|<Xt)> = |01) - ; ^|10). (25) 

y/2{e~« + 1) y/2(e-« + iy 

From Eqs. (25) we can see that with the time evolution p(t) tends to ||</?)(</?| + ||0}(<^| + -||00)(00|, which is a 
mixture of two-photon, one-photon and zero-photon zero-mode solutions, as analyzed in section (IV. A). Since the 
quantification of entanglement of mixed entangled qutrit states is not yet well understood, we can give an upper 
bound of the entanglement by the converxity property of entropy function, namely, 

E(p(t)) < {e ~* + l)2 E{\ V {t)){ip{t)\) + l ~ C 2 - E{\4>{t))W)\) + ^ E(\00)(00\) = E*{p{t)). 

In the above expression, E is defined by E = h( where h(x) = —xlog 2 x — (1 — x) log 2 (l — x) and C is the 

concurrence of the related density matrix. From Eq. (24) the two-photon pure state entanglement E(\ip(t)} (ip(t)\) 
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can be obtained. E(\<f>(t))((f)(t)\) is just a pure qubit entanglement, which can be calculated in the traditional way. 
Fig. 7 shows the time evolution of E*(p(t)) (solid line). One notes that the steady entanglement can be generated 
from the product initial state. 

The similar discussions can also be applied to the cases ofc=l,a = 6 = and b = l,a = c = 0. While the 
evolution of the upper bound E*(p(t)) for b = 1, a = c = is the same as the above, a slightly different evolution 
in the case ofc=l,a = 6 = 0is found, as shown in Fig. 7 by the dashed line. Nevertheless, the conclusion of the 
generation of the entanglement is same. 

Case 2, a = b=\,c= -1/V2 

In this case the solution is time-independent, 

p{t) = \y){ V \, 

where \ip) is the same as the form in section (IV. A). Once again we verify the existence of decoherence-free state. 
Case 3, a = b=\,c = 1/V2 

From the analysis of last subsection we know that environment can produce and destroy entanglement, which 
is dependent on the initial state. This is still true in the two-photon process. The Case 1 corresponds to the 
entanglement production. In this case we shall see the environment can destroy entanglement in certain initial state. 
The solution read 

p(t) = e-^\x)( X \+2( e -^^e-^)\v)(v\ + (l-e-^) 2 \00)(00\, 
|X) = i|02) + i=|ll) + i|20), 

H = >i> + >°>- 

We can see p(t) tends to zero-photon zero- mode solution |00)(00|. As a result, the entanglement is destroyed com- 
pletely. The time-evolution for E*(p(t)) is plotted in Fig. 7, as shown as dotted line. 

The discussion above is limited to very special initial states, which are formulated in the Fock states representation. 
It is very interesting to expand the present discussion to the Gaussian states, which can be created relatively easily and 
can be widely applied to quantum cryptography and quantum teleportation. In addition, the entanglement criteria 
[28-30] and the formation of entanglement [31] are also available. However, instead of the Fock state representation, 
it is very convenient to use the continuous variable representation to discuss the Gaussian states, which is beyond the 
present work. 



VI. SUMMARY 



In summary, we have investigated the master equation of two identical single-mode cavity fields coupled to a 
common quantum environment. Using the algebraic dynamical method, the two kinds of dynamical symmetries, 
namely the sp(4) D su(l, 1) D {Kq} and sp(4) D u(2) D {U, So} dynamical symmetries of the master equation have 
been explored. The first dynamical symmetry corresponds to the collective quasi-photon mode, while the second one 
is the dynamical symmetry for the real photons of the two cavities. According to algebraic dynamics, the system is 
intcgrable and its analytical solutions have been obtained. It is found that the system is affected by the environment 
through the rate operator consisting merely of collective quasi-photon operators and having collective quasi-photon 
eigen modes which are highly degenerate in terms of the real photon states. Based on the analytical solutions, the 
decoherence and entanglement properties of the quasi-photon zero mode of the rate operator have been studied in 
detail, and the surprising finding is that different entangled states can be produced from the disentangled initial 
states by environment effects in some cases. The results obtained are useful for the encoding scheme in quantum 
computation and for the entanglement generation of the 2-qubit and 2-qutrit cavity systems. 
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FIG. 1: The time evolution of the elements of density matrix with the initial condition of 




FIG. 2: The time evolution of the concurrence C(p) with the initial condition of a — 




FIG. 4: The time evolution of the concurrence C(p) with the initial condition of a = l/y/2. 
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FIG. 7: The upper bound E* of the entanglement of two-photon process for different cases: o = 1, b = c = (solid line), 
c — l,a = b — (dashed line), and a = b = |,c = 1/V% (dotted line). 



